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1. Introduction 

Present-day laser systems may deliver intensities up to 10 22 W cm -2 pQ at their focal 
spot. Even higher intensities of the order of 10 24 — 10 26 Wcm -2 are envisaged at 
the Extreme Light Infrastructure (ELI). Theoretical studies [2] suggested that in 
the interaction of a laser pulse with a thin foil, Radiation Pressure Acceleration 
(RPA) becomes the dominant mechanism of ion acceleration at intensities exceeding 
10 23 Wcm -2 . The radiation pressure dominated regime is attractive because of the 
foreseen high efficiency and because of the quasi-mono energetic features expected in 
the ion energy spectrum. Moreover, recent simulations suggest that multi-dimensional 
effects may allow a further increase of the ion energy [3]. 

At these extreme optical laser intensities / > 10 23 Wcm~ 2 , electrons become ultra- 
relativistic within a fraction of the wave period experiencing super-strong accelerations 
and therefore emitting relatively large amounts of electromagnetic radiation. Radiation 
reaction (RR) is the influence of the electromagnetic field emitted by each electron on the 
motion of the electron itself [4] and may become essential under the extreme conditions 
mentioned above. Early particle-in-cell (PIC) simulations [5] showed that RR effects 
become important at intensities exceeding 5 x 10 22 W cm -2 and increase nonlinearly 
with the laser intensity. 

In order to take RR effects self-consistently into account one should, in principle, 
solve the so-called Lorentz-Abraham-Dirac (LAD) equation [I] . It is well known that this 
equation is plagued by inconsistencies such as, for example, the appearance of "runaway" 
solutions in which an electron acquires an exponentially diverging acceleration even 
without any external field. However, it has been shown that in the realm of classical 
electrodynamics, i.e. neglecting quantum effects, the LAD equation can be consistently 
approximated by the so-called Landau-Lifshitz (LL) equation which is free from the 
mentioned inconsistencies [HE]. 

In this paper, we investigate RR effects in the interaction of a super-intense 
laser pulse with a thin foil in the RPA-dominant or "laser-piston" [2] regime by 
one- dimensional (ID) PIC simulations both for linear and circular polarization. Our 
approach is based on the LL equation of motion. We identify leading terms in 
the LL equation and discuss suitable approximations. On this basis we develop a 
straightforward numerical implementation of the RR force in a standard PIC code. 
PIC simulations with RR effects included have been previously performed for various 
laser-plasma interaction regimes by several groups, either using an approach similar to 
the LL equation [5] or using a different RR modeling [7HT2]. 

In our simulations, we check the RR's ability to reduce the electron heating which 
is responsible of the broadening of both the electron and ion spectrum. Indeed, recent 
studies for thick targets in the hole boring regime [H][l2] and ultrathin plasma slabs [9] 
suggested that the RR force cools the electrons and may improve the quality of the 
accelerated ion bunches. We found that in the linear polarization (LP) case, the peak in 
the energy spectrum has both a lower energy and a lower width when RR is included. At 
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the same time, the fraction of low energy ions is reduced. However, strong modulations 
appear in the ion energy spectrum after the acceleration phase both with and without 
RR and eventually the quasi- monoenergetic features are lost. In the circular polarization 
(CP) case, RR does not affect the ion energy spectrum significantly even at intensities 
of the order of 10 24 Wcm~ 2 . The differences between LP and CP appear to be related 
to the longitudinal electron oscillations driven by the J x B force in the LP case. These 
oscillations allow a deeper penetration of the laser pulse into the foil enhancing the 
effect of the RR force on electrons. In the CP case, significant RR effects are found only 
for laser and target parameters such that the laser pulse breaks through the foil due to 
nonlinear transparency, similarly to what was found in previous studies [9]. 

2. The Radiation Reaction force 

In classical electrodynamics, the effect of RR on the motion of an electron can be taken 
included by means of an additional force besides the Lorentz force. The additional RR 
force basically describes the loss of energy and momentum by an accelerated electron 
which radiates EM waves, so that the electron trajectory changes with respect to that 
predicted by the Lorentz force alone. In the LL approach [I] the RR force is written in 
a manifestly covariant form as 

f 1 = (d a F^u u u a ) + (F^F va u a + (F»%F ua u a ) U v) (1) 

6mc z Sm^c^ 

where m and e are the electron mass and charge respectively, w M = (7, yv/c) is its four- 
velocity and F^ v is the electromagnetic tensor relative to the total electromagnetic field 
acting on the electron except the field generated by the electron itself. 

The importance of RR effects on the electron motion depends on the strength and 
geometry of the EM fields, as well as on the electron energy which is generally a function 
of the amplitude and of the frequency of the field itself. One would thus need to know at 
least the scaling of the electron energy with the laser pulse parameters for a preliminary 
evaluation of RR effects as well as for a discussion on the limits of validity of the chosen 
theoretical approach and on suitable approximations to it. In the following discussion, 
we mostly refer to the case of the electron motion in a plane wave. For this problem, the 
LL equation has an exact analytical solution for arbitrary pulse shape and polarization 
of the plane wave [13] . Such solution thus provides a useful benchmark and reference for 
RR effects in superstrong laser fields. In a many-particle system such as a high-density 
plasma, the collective EM fields are generally much more complicated but the plane 
wave results may provide some guidance for their interpretation. 

We first recall that the LL approach is classical and quantum electrodynamics 
effects are neglected. In the interaction between an intense laser field (with peak 
intensity / and wavelength A) and an ultra- relativistic electron (with Lorentz factor 
of the order of 7) this is in general allowed if 7^/ /I cr "C 1 and 7A C /A <C 1 [4J, where 
I cr = cE 2 r j%n m 2.3 x 10 29 Wcm -2 is the intensity corresponding to the critical field 
E cr = m 2 c 3 /h\e\ of quantum electrodynamics [H] and A c = h/mc ~ 3.9 x 10 _7 /im is 
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the Compton wavelength. These conditions ensure that the momentum of the photons 
emitted or absorbed by the electron is negligible. Moreover, the force related to the 
electron spin might not be negligible in comparison to the RR force. In fact, the 
dynamics of a particle with a spin degree of freedom in an external electromagnetic 
field can be described in the classical framework by the Frenkel force [15] (see also [16J 
for a different derivation of both the RR and the spin force) 



ft' 

Js 



(2) 



where Q 1& 



m 



is the magnetic dipole moment four-vector and 



the Levi-Civita symbol (e 0123 = +1). The analysis of the case of a plane wave (electric 
field amplitude E, central frequency u and pulse length r) shows that the spin force is 
about ~ 7/a ~ 1377 times the term in the LL force (CQ) containing the derivatives of 
the field tensor, i.e. the term proportional to d\F^ v (here a = e 2 /he ~ 1/137 is the 
fine-structure constant). However, it can also be shown that the spin effects remain 
smaller than those due to the last term in Eq.([T]) if aa ur > 1 where a = \e\E/mujc 
(the effect of the last RR term cumulates with time). Since r > 2-k/uj and ao > 300 
in our simulations the latter condition is well satisfied. It is therefore consistent, in a 
regime where RR effects are relevant and quantum effects are subdominant, to neglect 
both the spin force and the first term of the RR force in Eq.(TjQ). 

The PIC simulations with the RR force included are performed in the laboratory 
frame, i.e. the frame where the plasma target is initially at rest. In the laboratory frame 
we write down the LL equation in three-dimensional, non-manifestly covariant form as 
dp 



dt 
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(3) 



where p is the electron momentum, r e = e 2 /mc 2 ~ 2.8 x 10 -9 fim is the classical electron 
radius, A = 2ttc/u is the laser wavelength and we use dimensionless quantities as in the 
PIC code. Time is normalized in units of uj~ 1 , space in units of cw _1 , momenta in units 
of mc. Consequently, EM fields are normalized in units of muc/\e\ and densities in 
units of the critical density n c = muj 2 / Aire 2 . The first RR term of Eq.([n]), i.e. the one 
containing the "total" time derivative of the EM fields, corresponds to the negligible 
term in the manifestly covariant LL Eq.flTJ and is reported here for completeness but 
neglected in the calculations for the above explained reasons. 

Since RR effects are important for ultra-relativistic electrons 7 ^> 1, the last term 
in Eq.(j3]) (proportional to 7 2 ) dominates over the preceding one. From a practical point 
of view, the smaller term may often be neglected even though the on-shell condition 
1 is lost neglecting this term. Although single particle and PIC tests with and 
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without this term showed no significant difference, both terms were included in our PIC 
simulations. It is possibly instructive however to neglect for a moment the smaller term 
and write down an effective reduced LL equation in the lab frame 

£ = fe -<v (4, 
where fj, = — (E + v x B) and d is given by 



4 r e 



3 A 

4 r e 

-7T— 

3 A 



d = ( -vr^ )<f E 2 — (v • E) 2 + v 2 B 2 - (v B) 2 - 2v ■ (E x B 



7 2 [f 2 - (v ■ f L ) 2 ] > 0. (5) 



In Eq.flH]), RR effects appear as a "friction" term with a nonlinear and anisotropic friction 
coefficient given by d. When an electron "crosses" an EM field, it feels a viscous force 
opposite to its velocity. 

For an ultrarelativistic electron, the friction coefficient d may be used as a measure 
of the strength of the RR force in units of muc. In the case of motion in a plane wave, 
d may be compared directly to the normalized wave amplitude ao- Setting E x B along 
the positive x axis, the RR force vanishes (d — > 0) when v x — > 1, has its maximum value 
(d ->■ (|vr^) 7 2 4ao) when v x ->■ -1 and finally d -)■ (§7r^-) 7 2 Oq when (w 2 + v 2 ) -)■ 1. 

The friction effect of the RR physically corresponds to the incoherent emission of 
high frequency radiation by ultrarelativistic electrons. When the RR is included in the 
numerical simulation of a collisionless, relativistic plasma, it is typically not feasible to 
resolve electromagnetic waves at such high frequencies, much larger than the inverse 
of the temporal resolution. Thus, it is assumed that such radiation escapes from the 
system without re-interacting with other electrons. Notice that even a solid-density 
plasma is transparent to such radiation, since in the present regime the RR effect is 
mostly due to the emission of radiation with photon energies in the MeV range, while 
the plasma frequency corresponds to at most a few hundreds of eV. From the point 
of view of energy balance, then, the energy radiated at high frequencies appears as a 
loss term or "dissipation" . The percentage of radiative loss is measured by comparing 
the energy balance simulations including RR with simulations without RR, where the 
total energy of fields and particles is conserved within the limits of numerical accuracy 
(typically within 1% in our PIC code). 

It may be worth recalling that, for what concerns the LL equation of motion, 
energy and momentum are not conserved exactly for the single electron. This is due to 
some terms which are neglected when deriving the LL equation from the LAD equation 
under the assumption that the radiation force in the instantaneous rest frame of the 
electron is much smaller than the Lorentz force [3]. However, the neglected terms are 
much smaller than quantum corrections jl], thus the approximation is consistent with 
a classical treatment. In a different approach to the RR force recently presented in Ref. 
[IT] , a different couple of classical equations of motion are derived phenomenologically 
starting from the requirement of energy-momentum conservation of the system of the 
electromagnetic field plus the radiating electron. 
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3. The numerical approach 

Our PIC code is based on the standard, widely used Boris particle pusher [18] and 
leap-frog schemes to advance and accelerate particles. We developed a simple numerical 
scheme to insert the RR force in the PIC code while keeping the standard particle pusher 
for the Lorentz force unchanged. As will be clear below, this scheme is based on the 
assumption that the acceleration of particles is dominated by the Lorentz force, with 
the RR force giving a smaller, albeit non negligible contribution. 

We write the total force f acting on the electron as the sum of two forces (already 
introduced) and with 

f R = - Qvr^ {f L x B - (v ■ E)E + 7 2 [f 2 L - (v ■ E) 2 ] v}. (6) 

Then, the equation of motion of the electron reads 

f-f-* + * (T) 

Assuming that forces and momenta are known at integer and half-integer timesteps 
respectively, the full leap-frog step is 

D (n+l/2) _ D (n-l/2) 

^ = f (n) = t ] + # ) (8) 

where At is the timestep. Now, we consider the leap-frog step for two "helper" momenta 
p L and p R 

(n+l/2) (n-1/2) _(»+l/2) _(n-l/2) 

Pl ~ Pl _ f («) Pr ~ Pr = f (n) , Q s 
At L ' At R 1 ' 

and assume = 1 ' 2 ^ = p( n_1 / 2 ). Thus, from the above equations we easily 

obtain 

p (n + l/2) = p (n + l/2) + p (n+l/2) _ p (n-l/2) (1Q) 

This means that, starting at time t^> and position with p( n_1 / 2 ) ; firstly p(™ +1 / 2 ) 
and p^ t+1 ^ 2 ' ) are calculated independently using and respectively, and finally 
Eq. ffTU]) is employed to obtain the full leap-frog step p( n+1 / 2 ). It is worthwhile noticing 
that this is a general result as we have used only the superposition property of the force 
without any assumption about fi, and f^. 

The previous algorithm allows to keep the standard leap-frog pusher for the Lorentz 
force and to develop an independent pusher for the RR force alone. Using Eq. (Q we 
can recast Eq. ( JTUl) as 

p (n + l/2) = p (n+l/2) + f (n) At = p (n + l/2) + ^ (n) 

Now, in order to compute the momentum change from step n — 1/2 to n + 1/2 due 
to the Lorentz and the RR force, an estimate of the electron's velocity at halfstep n is 
needed. To this aim, we first advance p( n-1 / 2 ) to p^ +1 / 2 ) using the Boris pusher for the 
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Lorentz force, then we use p^ 
v( ra ) at half time step as 



(n+l/2) 



to estimate the total momentum p*™) and velocity 
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1 + (pW) 2 
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Next we use Eqs. (fT2l and (fT3"|) together with the fields E^ n \ B^ n ^ at half time-step to 
compute the full term fjj according to Eq. ([6]). This task is particularly simple because 
many terms of f# can be written by directly (see Eq. fl6])). 

This particle pusher was tested comparing the numerical results for a single electron 
in a monochromatic plane wave both with the known analytical solution [13] and with 
the numerical solution obtained using a 4th order Runge-Kutta scheme. These numerical 
calculations confirmed that the inclusion of the RR force according to the above method 
preserves the accuracy and stability of the standard Boris pusher algorithm. The range 
of intensities in the tests was from 10 22 W cm -2 to 10 24 W cm -2 with A = 0.8 fim. Taking 
as an example case an electron with initial momentum p xo = —200 mc and a wave with 
ao = 350 and A = 0.8 /zm, we found our particle pusher to yield a phase error in the 
longitudinal momentum of ~ 0.1(2tt/oj) after a run time of 500 u~ l using a timestep 
At = 0.01 uj~ 1 . The corresponding relative error in the displacement in the direction 
of wave propagation was ~4x 10~ 4 . The one-particle tests were performed using the 
complete expression of the LL force ([3]) with the fields and their derivatives as given 
functions of space and time. These tests also confirmed that the derivative term in the 
LL force ([3]) is negligible. The inclusion of the RR force in the PIC code according to 
the above described approach leads to approximately a 10% increment of the computing 
time. 

4. The PIC simulations 

We performed PIC simulations with a plasma slab of ions (protons) with uniform initial 
density uq. Since our primary aim is to evaluate the importance of RR effects on laser- 
plasma dynamics and ion acceleration in the regime of radiation pressure dominance, 
we restrict ourselves to a one-dimensional (ID) geometry for the sake of simplicity and 
the possibility of using high numerical resolution. Mult i- dimensional effects, which 
may be important to determine the features of ion acceleration in this regime [3], 
will be presented in forthcoming publications; preliminary 2D simulations [19] show 
qualitatively similar trends to the ID case. It is worth noticing that, as the momentum 
space in the ID PIC code is already three-dimensional, our numerical approach can be 
readily implemented in a multi-D code employing the same particle pusher. The modest 
increase in computational time implied by our method might be essential to be able to 
perform large-scale multi-D simulations with RR included. 
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We report results for a laser intensity / = 2.33 x 10 23 W cm -2 and a laser wavelength 
A = 0.8 fim, corresponding to a dimensionless parameter oq = 328. In all the simulations, 
the density uq = 100n c and the profile of the laser field amplitude has a "trapezoidal" 
shape in time with one-cycle, sin 2 -function rise and fall and a five cycles constant plateau. 
The laser pulse front reaches the edge of the plasma foil at t — 0. The foil thickness 
is i = 1A in all the simulations except for the "transparency" case reported below, for 
which i = 0.3A. We considered both Circular (CP) and Linear (LP) polarization of the 
laser pulse. The parameters are similar to those of the 3D simulations in Ref. [2] where 
the laser pulse was linearly polarized. According to [2] , RPA dominates the acceleration 
of ions in the plasma foil when the laser intensity / > 10 23 Wcm~ 2 . To our knowledge, 
neither RR effects nor CP have been studied so far in such a regime of laser and plasma 
parameters. The effects of CP have been studied extensively at lower intensities (see 
e.g. Ref. (20] and references therein) showing that, with respect to LP, the use of CP 
quenches the generation of highly relativistic electrons making RPA dominant also at 
such lower intensities. Concerning RR effects, in Ref. [2] it was suggested that the higher 
velocity the plasma foil is accelerated to, the lower the RR force becomes because of the 
relativistic increase of the laser wavelength A' in the foil frame, making the RR strength 
parameter ~ r e /\' increasingly small. The expected quenching of RR effects may also 
be explained with the help of the "reduced" LL equations d!|)-(j5]): when the foil moves 
coherently with a velocity close to c, the amplitude of the reflected wave is strongly 
reduced at any time in the laboratory frame; thus, the electrons at the surface of the 
foil can be considered as moving with a velocity v x ~ c in the field of the incident plane 
wave and parallel to its propagation direction, and the RR force almost vanishes. 




Figure 1. Ion energy spectrum f{E) at t — 46 T with (red) and without (black) 
RR for CP. The laser intensity is / = 2.33 x 10 23 Wcm" 2 and the target thickness is 
i = 1A. See the text for the parameters common to all the simulations. 



In the CP case, we found that RR effects on the ion spectrum (distribution of 
protons per unit energy) are negligible as shown in Fig. [1] for a time t = 46T where 
T = A/c is the laser period. Even at higher intensities, RR effects on the ion spectrum 
are weak provided that there is not a strong transmission of the laser pulse through 
the foil. In the simulation corresponding to Fig. (U the laser pulse penetrates into the 
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Figure 2. Snapshot at t = 1.7 (T) of the "skin"layer of the foil for CP and 
J = 2.33 xlO 23 Wcm" 2 . The foil was initially placed between a; = 14 A and a; = 15 A. a) 



The electron density (black), the modulus of the transverse electric E±_ = J E 2 + E% 
(blue) and magnetic B^ = « + B| (red) fields. Distribution of b) longitudinal 



momentum p x , c) modulus of the transverse momentum pj 
coefficient d. 



- p 2 z and d) friction 



plasma for a small distance of the order of A/20, and the fields in the plasma are much 
smaller than the fields in vacuum. As a consequence, the friction coefficient d introduced 
in Eq. (JSJ) is very small compared to ao- The spatial profiles of both the fields and the 
coefficient d in the "skin" layer are shown in Fig. |2] The order of magnitude of the 
normalized transverse momentum is ~ 10 and of the friction coefficient is d ~ 10~ 2 . 
It is worth mentioning that Fig. |2] shows a snapshot at t = 1.7 (T) but the typical 
values of the friction coefficient d are always of the same order of magnitude for CP. 
In contrast, for LP the friction coefficient d attains much larger values at the same 
instant, as discussed below. We also notice that, for CP, we obtain qualitatively similar 
results also at higher intensities, up to 10 24 Wcm~ 2 . However, at such extremely high 
intensities the condition of validity of the classical approach {^f\J E / E CT < 1) tends to 
be violated, so at least such results should be taken with caution and an analysis based 
on quantum RR effects might be necessary. 

Reducing the foil density or thickness, the laser pulse may break through the foil. In 
this case more electrons move in a strong electromagnetic field becoming ultrarelativistic 
in a fraction of wavecycle and RR effects strongly affect the ion spectrum, as shown in 
Fig. |3j In particular, when RR is included, peaks in the energy spectrum appear at 
energies higher than in the case without RR. This result is similar to that obtained in 
Ref. [9] at lower intensities (~ 10 22 Wcm~ 2 ), where it was suggested that RR effects 
"improve" the ion spectrum in the optical transparency regime. Our explanation is that 
the effective "dissipation" due to RR leads to a later breakthrough of the laser pulse 
through the foil, favoring a longer and more efficient RPA stage. However, comparing 
Fig. |3]with the thicker target case in Fig. [lj it is evident that the spectrum becomes very 
far from monoenergetic, while the maximum ion energy increases only slightly. Hence 
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in our simulations "optimal" conditions for ion acceleration are found for the case of 
Fig. [TJ for the corresponding laser and plasma parameters, RR effects are negligible. 




0.0 0.5 1.0 1.5 2.0 

E (GeV) 



Figure 3. Ion energy spectrum f(E) at t = 46 (T) for a simulation with the same 
parameters as Fig. [T]but with a target thickness £ = 0.3A. In this case the laser pulse 
breaks through the foil and RR effects are evident. 

In the LP case the foil is accelerated by radiation pressure too but, unlike the CP 
case, the laser pulse does penetrate up to a fraction of the order of A/4 at the front 
surface of the foil, as shown in Fig. HI The two snapshots are selected both to show 
values of d close to its maximum in time and to make a direct comparison with the CP 
case of Fig. El It is found that a larger fraction of electrons at the front surface move in 
a strong electromagnetic field of the same order of the vacuum fields. In this case, the 
friction coefficient function d reaches values of d ~ 10 2 (Fig. Hj) which are comparable 
with the Lorentz force (a = 328). The deeper penetration of the laser pulse is correlated 
with the strong longitudinal oscillatory motion driven by the oscillating component of 
the J x B force which is suppressed for CP. Large numbers of electrons are pushed 
periodically inside the foil producing strong fluctuations of the electron density (see 
Fig. H] part a)). 

For LP, the ion energy spectrum is significantly affected by RR effects. The 
spectrum is fairly peaked with a smaller energy spread and lower peak energy than in 
the case without RR (Fig. [5]). In general, as observed in many simulations the spectral 
peak produced by RPA broadens with increasing electron "temperature", since hot 
electrons drive the expansion of the plasma leading to additional, non-monoenergetic 
ion acceleration. The smaller energy spread observed when RR is included can be 
thus traced back to the radiative cooling of the most energetic electrons. Moreover, a 
significant fraction of ions on the low energy tail of the spectrum is observed without 
RR, but disappears when RR is included. The fractional difference in the ion energy 
with vs without RR is of the order of the fraction of the laser pulse energy that is "lost" 
as incoherent emission (Fig. [6]). For / = 2.33 x 10 23 Wcm -2 , about 20% of the total 
pulse energy is lost as incoherent radiation (Fig. [6]). 

When RR is "switched off", part of the "skin" layer of the foil is left behind and 
a significant fraction of ions is present on the low energy tail of the ion spectrum 
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Figure 4. Snapshot at t = 1.6 (T) and f = 1.7 (T) of the "skin"layer of the foil for 
LP and I = 2.33 x 10 23 Wcm -2 . The strong longitudinal oscillations driven by the 
J x B force allow a deeper penetration of the laser pulse into the foil compared to the 
CP case. The foil was initially placed between x = 14 A and x = 15 A. a) The electron 
density (black), the modulus of the transverse electric \E y \ (blue) and magnetic \B Z \ 
(red) fields. Distribution of the b) longitudinal momentum p x , c) modulus of the 
transverse momentum p± = \p y | and d) friction coefficient d. Notice the change of the 
scale from t = 1.6 (T) and t — 1.7 (T) in the frames b), c) and d). We remark that the 
longitudinal momentum p x distribution changes of orders of magnitude in 0.1 (T) due 
to the J x B force. 
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Figure 5. Ion energy spectrum at t — 12 (T) with (red) and without (black) RR for 
LP and / = 2.33 x 10 23 WcuT 2 . 



(Fig. [5]). To explain this effect, we first recall that in the first stage of RPA two ion 
populations may be produced, corresponding to a coherently moving "sail" and to a 
trailing "tail" [20]. Ions in the tail will eventually remain behind the sail if their charge 
is neutralized by returning electrons; otherwise, they will be accelerated by their own 
space-charge field and may move to the higher energy side. When the foil is still non- 
relativistic in the laboratory frame, the RR force has larger values when the electrons 
counter-propagate with respect to the laser pulse and therefore the electron backward 
motion is strongly impeded when the RR force is included. This effect prevents an 
efficient neutralization of the ion charge in the tail by returning electrons, explaining 
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Figure 6. Fractional energy absorption in function of time for LP and intensity 
I = 2.33 x 10 23 Wcm~ 2 . Electron kinetic energy with RR (blue) and without RR 
(black), ion kinetic energy with RR (orange) and without RR (red) and the fraction 
of energy lost in system (green). 
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Figure 7. Ion energy spectrum at t = 46 (T) with (red) and without (black) RR for 
LP and / = 2.33 x 10 23 WcuT 2 . 



Equation (j2J) suggests that the RR force is mainly a nonlinear friction force. For 
/ = 2.33 x 10 23 Wcm~ 2 , about 20% of the total pulse energy is "dissipated" by the 
RR force during the laser-foil interaction (Fig. [6]) which lasts about 22 cycles (30 cycles 
without RR). As stated previously, such "dissipated" energy accounts for the incoherent 
radiation escaping from the plasma. During the laser-foil interaction, such flux of 
incoherent radiation shows itself in a missing pulse energy while ions have almost the 
same total energy in both cases and their spectrum is quasi- monochromatic (Fig. |5]). 
However, after the acceleration phase by the radiation pressure of the laser pulse, a 
20% of missing pulse energy implies about the same amount of missing final ion energy 
(Fig. |6]). Moreover, a significant fraction of hot electrons are produced by the J x B 
force. Such electrons can drive an expansion of the foil, strongly increasing the ion 
energy spread after the laser- foil interaction phase (Fig. [7j). 
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We remark that just changing the laser polarization from CP to LP, the friction 
coefficient d increases by up to four-orders of magnitude due the enhanced laser pulse 
penetration in the foil by the J x B-driven longitudinal oscillations. Then, the electrons 
move in a strong electromagnetic field becoming ultrarelativistic and the "friction" term 
of the RR force becomes non negligible. These results are a relevant test of the conjecture 
in Ref. j2] that RR effects would be weak as the foil motion becomes relativistic. Our 
simulations suggest that this picture strictly holds only in the CP case, where almost 
all of the foil moves at relativistic speed in the same direction of the laser pulse. In 
the LP substantial fraction of electrons has both an ultrarelativistic motion in 

the transverse direction and a strong oscillatory motion in the longitudinal direction, 
leading to significant RR effects. 

The dependence of RR effects on the pulse polarization was also studied in Refs. 
[T014T2] for thick targets ("hole boring" regime of RPA) and long pulse durations. It 
was also found that RR effects are stronger for LP, although they are not negligible for 
CP [12J. These results cannot be compared straightforwardly to our findings because 
of the quite different laser and plasma parameters, leading to a different dynamics. For 
instance, in the thick target case the laser-plasma surface oscillates also for CP ( "piston 
oscillations" [12]) and a return current of electrons counterpropagating with respect to 
the laser pulse is generated; this effect is likely to enhance radiative losses. 

5. Conclusions 

We summarize our work as follows. Radiation Reaction effects on Radiation Pressure 
Acceleration of plasma slabs by ultraintense laser pulses were studied by one-dimensional 
PIC simulations. The RR force was included via the Landau-Lifshitz approach. The 
numerical implementation allows the addition of RR effects to any PIC code based on 
the standard Boris pusher algorithm for the acceleration of the particles, at a small 
computational cost. 

We compared results for Circular and Linear Polarization of the laser pulse. For 
CP, we found that RR effects become relevant only for plasma targets thin enough to let 
the laser pulse break through the foil. In this case the inclusion of RR effects leads to an 
increase of the ion energy. Such increase is however not very significant with respect to 
a case with the same laser parameters but a thicker target, for which the breakthrough 
of the laser pulse does not occur and RR effects are negligible. 

For Linear Polarization, we found that RR effects are significant, leading to some 
tens of percent of energy loss by incoherent emission and to a reduction of the peak ion 
energy by a similar percentage. Although RR effects produce a somewhat more peaked 
energy spectrum during the acceleration stage, the final spectrum is anyway dominated 
by a post-acceleration evolution, presumably driven by high energy electrons. 
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